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Abstract 



The hyperfine interactions of the constituent quark model provide a nat- 
ural explanation for many nucleon properties, including the A — N splitting, 
the charge radius of the neutron, and the observation that the proton's quark 
distribution function ratio d{x)/u{x) — ► as x — ► 1. The hyperfine-perturbed 
quark model also makes predictions for the nucleon spin-dependent distribu- 
tion functions. Precision measurements of the resulting asymmetries A\(x) 
and A™(x) in the valence region can test this model and thereby the hypothesis 
that the valence quark spin distributions are "normal". 
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I. INTRODUCTION 



The quark model has enjoyed so much success as a qualitative guide to hadronic structure 
that the discovery that only about 30% of the proton's spin could be attributed to quark 
spin came as a surprise. Since the quark model remains unjustified within QCD, it is a 
misnomer to call this "proton spin surprise" the "proton spin crisis" . However, whatever we 
call it, this result has generated much very productive experimental and theoretical activity. 

While in general the spin of the proton could reside on any mixture of its quark and 
gluon constituents or in their orbital angular momenta, a conservative interpretation [|TJ of 
the current situation is that the valence quarks carry the spin expected by the quark model 
but that the low x sea of qq pairs is negatively polarized. In this case £ (defined to be twice 
the expectation value of the quark plus antiquark spin along the spin direction of a polarized 
proton, so that £ = 1 would saturate the proton spin), when decomposed into its valence 
and sea components, would be 

S = S l) + ^A(g + g) sea (1) 
q 

where £„ = J dxT, v (x) is twice the spin on the valence quarks and Aq sea = J dxAq sea (x) 
and Aq sea = J dxAq sea (x) are, respectively, twice the spin on the sea quarks and antiquarks 
of flavor q. If the valence quarks were in nonrelativistic S'-waves as in the naive quark 
model, then S t , would be unity. However, as has been appreciated for nearly thirty years 
||, in realistic valence quark models lower components of quarks spinors convert about 
25% of the quark spin into orbital angular momentum so that £„ ~ 0.75. If in addition 
each of the three light quark flavors carries A(q + q) sea ~ —0.15, a very modest per flavor 
effect, £ ~ 0.30 would follow. Sea quark polarizations of just this sign and magnitude 
have recently been obtained in a realistic model of qq pair creation [p]] . (In a more general 
context, such small A(q + q) sea values are perfectly consistent with a 1/N C expansion of QCD 
(where sea quarks appear at order 1/N C via quark-antiquark loops). Note that the condition 
A(q + q) sea <S 1, not how accurately approximates unity, determines the applicability 
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of the 1/N C expansion: any nonzero A(q + q) se a would lead to a "spin crisis" as Nf (the 
number of light flavors) tends to infinity. ) 

In the conservative scenario just described, both the 25% relativistic quenching of spin 
from E„ and the negative polarization of A(q + q) se a are compensated by orbital angular 
momentum. In general, however, we are only guaranteed that 

£„ + A (<? + Q)sea + 2L q + = 1 (2) 

(where L q is the quark and antiquark orbital angular momentum and |E 9 is the total angu- 
lar momentum residing in the gluonic fields), so major experimental efforts are planned to 
measure the component parts of Eq. (2) in an effort to disentangle the "spin crisis" . These 
efforts begin with planned extensions of deep inelastic lepton scattering measurements of 
the proton and neutron spin structure functions down to very small x to complete the inte- 
grals required to calculate E, and studies of the Q 2 -dependence of spin structure functions 
to make inferences about Ag(x), the gluon helicity contribution to S 9 (x). Major efforts 
are also planned to directly measure Ag(x) based on helicity-dependent gluon-parton cross 
sections. In addition to these classical inclusive measurements, flavor-tagging semi-inclusive 
experiments are planned to measure separately As sea (x), As sea (x), Au sea (x), Ad sea (x), and 
also the quark contributions Au(x) = Au v (x) + Au sea (x) and Ad(x) = Ad v (x) + Ad sea (x). 
(Note that it is not possible to experimentally separate the quark contributions Au sea (x) 
and Ad sea (x) from Au v (x) and Ad v (x): this separation is conceptual only.) Additional 
complementary information on the ss content of the proton is expected from planned mea- 
surements of the electric and magnetic form factors G S E and G S M of the s7 M s current using 
parity-violating electron-nucleon elastic scattering. 

Given the substantial effort being devoted to this problem, it is surprising that we still do 
not know whether our original simple picture of the spin structure of the valence quarks is 
right! To some degree this is because this question is not well-defined: in contrast to other 
methods (e.g., QCD sum rules 0), the quark model is not normally embedded in a field- 
theoretic framework. As a result, there are many difficulties in making comparisons between 
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the "predictions" of the quark model and the precisely defined quantities measured in deep 
inelastic scattering. As two illustrations of such difficulties, I note that: 1) the separation 
of Eq. (2) is Q 2 -dependent (e.g., Ag might be small at low Q 2 but very important at large 
Q 2 ) and, as mentioned above, 2) the u and d contributions to A«„(x) and Ad v (x) cannot 
be disentangled from those to Au sea (x) and Ad sea (x). However, beyond x ~ 0.3, sea quarks 
and antiquarks are scarce and, since gluons are too, such intrinsically field-theoretic issues 
as the factorization scheme dependence of E„(x) associated with the gluon anomaly || may 
be neglected. Thus while the integral values Au v = J dxAu v (x) and Ad v = J dxAd v (x) 
cannot be checked, those fractions of the distributions Au(x) and Ad(x) extending beyond 
x ~ 0.3 may be compared with valence quark model expectations with a residual ambiguity 
associated only with their Q 2 evolution. Although in what follows I will imagine distribution 
functions devolved to the "quark model scale" Ql ~ 1 GeV 2 , given this residual ambiguity 

1 will avoid predictions of the x-dependence of distribution functions and focus instead on 
the polarization asymmetries A1(x) and A™(x) which depend only on ratios of distribution 
functions and which should therefore have minimal Q 2 -dependence. It is unfortunate that 
the current experimental situation for A\(x) and A™(x) for x > 0.3 leaves much to be desired 
(see Figs. 1): it is even consistent with the naive SU(6) predictions. 

What are the valence quark model predictions for the resulting polarization asymmetries 
A\(x) and A^x) in the valence region? Ignoring Q 2 evolution, they are 

_ AAu v (x) + Ad v (x) 
l[X) ~ 4u v (x) + d v (x) [6) 

= 4Ad„(a;) + Au v (x) ^ 
Ad v (x) + u v (x) 

where u v (x) and d v (x) are the unpolarized valence distribution functions which integrate to 

2 and 1, respectively. From these formulas it is clear that the predictions depend on knowing 
the interplay between the valence quark spin and momentum wavefunctions so that there 
can be no unique prediction of the valence quark model for these asymmetries. However, I 
will argue here that its predictions are sufficiently well-determined that they can be used to 



answer the simple question of whether the valence spin structure is "normal" or not. 

Aside from this observation, there is little in this paper that could not be extracted 
from earlier work on this subject to which I will refer below. However, the results of this 
earlier work vary widely since they are based on diverse methods of dealing with rela- 
tivistic internal quark motion, various prescriptions for boosting to the infinite momentum 
frame, ad hoc versus dynamical origins for the assumed S , f/(6)-breaking, potential versus 
bag models, and choices of quark masses. Here I will assume that the hyperfine interac- 
tion is responsible for the d(x)/u(x) ratio as x — > 1, and then normalize predictions for 
the valence quark spin distribution functions to the data on this ratio. In doing so, I 
will not only avoid much model dependence, but also most of the pitfalls discussed above 
associated with not knowing precisely how to embed the quark model in field theory. 
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Figure 1(a): Data [|J on A\ and the prediction (shaded band) of the model described in 
the text; the SU(6) prediction is | (dotted line). 
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Figure 1(b): Data || on A" and the prediction (shaded band) of the model described in the 
text; the SU(6) prediction is 0. 

The body of this paper builds up to these predictions in steps. In the next Section I 
will review the naive SU(6) predictions and then modify them within the context of SU (6) 
by allowing the quarks to have relativistic internal motions. I then describe the breaking 
of S'[/(6)-symmetric quark spin distributions in the hyperfme-perturbed quark model and 
close with a brief historical overview. 

II. THE SU(6) AND "RELATIVIZED S77(6)" DISTRIBUTION FUNCTIONS 

I begin by recalling that in SU(6) one may simply write 

p T= uudC A i) S x X + (5) 

n t= dduC A ip s xl (6) 

where since C A and ip s are the antisymmetric color and symmetric L = spatial wavefunc- 
tions, 
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x x + = -y^(ut + m -2 m) (7) 

is the unique spin-1 wavefunction which is symmetric in the first two quarks as required by 
the Pauli principle ||. In the nonrelativistic SU(6) quark model one therefore expects 
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u v T (x) = ~v S u(6)(x) (8) 

u v I 0) = g«si7(6)(aO (9) 

d v T (^) = g«si7(6)(a5) ( 10 ) 
2 

d « I 0*0 = g w 5(7(6) (ar) (11) 

where t's!7(6)( a; ) is the universal SU(6) distribution function associated with ip s . These 

distributions lead to the standard SU(6) predictions d(x)/u(x) = 1/2, A\(x) = 5/9 and 
Ai(x) = 0, and Ga = 5/3. When the relativistic quenching mentioned above @ is turned 
on, it creates an x-dependent probability which we denote by \ca{x) for a spin up (down) 
quark to be flipped to down (up). This reshuffling of probability leads to the "relativistic 
S77(6)" spin distributions 

5 2 

u v ] {x) = [-- -c A {x)}v su{%) (x) (12) 

1 2 

U V l (S) = [g + gC A (x)]^ (6) (x) (13) 

4 T (ac) = [- + -CA(a:)]usr/(6)(:c) (14) 

2 1 

^ I ( x ) = [g - gCx(^)]vflU(6)(ac) (15) 

where with the nucleon expectation value (1 — c^(a;))jv — §GU — 0.75, the integrated valence 
spins become 



Au v ~ +-G A 



(16) 




(17) 



so that the "relativistic S77(6)" spin distributions satisfy the Bjorken sum rule. However, 

among other problems, the "relativistic S77(6)" model still makes the incorrect prediction 
d(x)/u(x) = 1/2. Note that the model also predicts that A\{x) = [1 — ca(x)]5/9 and 
A™{x) — as x — > 1, which, since c^(rr) — > as as — > 1, is not obviously wrong (see Fig. 1). 

III. PREDICTIONS OF THE HYPERFINE-PERTURBED QUARK MODEL 

Since the zeroth-order nucleons are pure S-waves, in the hyperfme-perturbed quark model 
0, only the Fermi contact part of the hyperfine interaction (the Si-SjS 3 {fij) force responsible 
for the A — N mass splitting) is operative in perturbing the nucleon's energy in first order. 
What does this perturbation do? In the nucleon rest frame, quark pairs with spin 1 have 
their energies raised (as in the A) while pairs with spin zero have their energies lowered. 
Since x X nas the two u quarks in a pure spin one state, while each ud pair is in a mixture of 
spin one and spin zero (with spin zero dominant so that the net perturbation in a nucleon 
decreases its energy), up quarks acquire higher average energy than down quarks. This 
physics then immediately suggests that the neutron will have a negative charge radius and 
that d(x)/u(x) will vanish as x —>■ 1 [jSHT0|1 . Since the individual spin components of x X are 
not in an eigenstate of the hyperfine interaction, it is less obvious what the effects are on 
the spin-dependent distribution functions. 

These effects are encoded in the L = component of the hyperfme-perturbed wavefunc- 
tion 



where (ip p ,ip x ) are mixed symmetry wavefunctions of the permutation group S3 which are 



uudC A cos6 m ^ s X X + + sin6 m J -{Vx+ ~ ^ X xX) 




(18) 
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antisymmetric (p) and symmetric (A) under 1 2 interchange, where 



:(U - IT) T 



(19) 



and where O m is a small mixing angle induced by S77(6)-breaking interactions |9|]. (Since, 

as explained above, the L = ground state energies are perturbed in first order only by the 
Si ■ Sj interaction, one can ignore L = 2 and totally antisymmetric L = admixtures.) It 
follows that the rest frame probabilities of spin up and spin down d quarks are, to first order 



in Or, 



p(dT) 



(20) 



P(di) 



0" 



(21) 



since the ip p x+ piece of the wavefunction does not interfere with the other terms in the 
probability distribution. In these formulas I have suppressed coordinate labels which indicate 
that the probability P(d f) (P(d |)) is that for finding a spin up (spin down) d quark at a 
point r<i while the two up quarks are at positions a and b. 
Similarly one finds 



p(«T) 



■0" 



(22) 



P(ui) 



2 2 



S n UP 



(23) 



where now the wavefunction ip p does play a role. I have now suppressed coordinate labels 
which indicate that the probability P(u f) (P(u I)) is that for finding a spin up (spin down) 
u quark at a point f u while the other up quark is at position a and the d quark is at position 

0. 

Note that, as advertized, the net leading-order effect of the >S[/(6)-breaking in the spin- 
averaged probabilities is to create distributions of mixed symmetry that allow the d quark 
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to have a different probability distribution from the two u quarks. With the calculated 
quark model value |J sin6 m ~ —0.23, the distortion of the S'[/(6)-symmetric probabilities 
is substantial. I now make the natural assumption that this distortion translates into the 
observation that d(x)/u(x) — ► as x — > 1, and associate the measured u(x) and d(x) with 



functions u v (x) and d v (x) associated with the spin-averaged probability 
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This remarkably simple picture then leads to the "standard" prediction Fg/F^ — > \ as 
x — > 1. 

The predictions of Eqs. (20)-(23) for A\(x) and A^ix) may easily be deduced using the 
properties of the mixed symmetry pair of wavefunctions (tp p , ip x ) under the permutation 
group S3. Such an analysis reveals that the hyperfine interactions have distorted the distri- 
butions of u I, d |, and d I identically, and that the entire dominance of u quarks as x — > 1 
is due to u f (x). This means that 

U v t (x) = [1 - ^C A (x)]u v (x) - -[1 - C4 (£)]<,(£) (24) 

m„ I (x) = i[l - c A {x)]d v (x) + ^c A (x)u v (x) (25) 
it(x) = -[l + -c4i)K(i) (26) 

4 I (x) = -[1 - -c A {x)\d v {x) (27) 

The resulting predictions for A\(x) and A"(x) in the valence region, shown in Fig. 1, can be 
obtained without engaging in an elaborate parameterization of structure functions. Using 
the rough parameterizations d(x)/u(x) ~ k(1 — x) as x — >• 1 (with 0.5 < k < 0.6) and 
ca( x ) — nx(l — x) n (which builds in c A (x) — > as x — » 1 and x — > and for 2 < n < 4 
gives the required quenching of G A ), produces the narrow bands shown in the figure. As 
x — > 1 both A\ and A" tend to 1, but I show the predictions only in the region where the 
valence quark wavefunction is large since very small effects might become important at the 
endpoint [jTT|| . 
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IV. SOME HISTORY 



The history of the prediction of the effects of S77(6)-breaking on the quark distribution 
functions in the valence region is somewhat convoluted. It perhaps begins with the parton 
model discussion by Feynman [[12] who argues that as a u or d quark approaches x = 
1, it must leave behind "wee" partons with either I — or I — 1, and that these two 
configurations are unlikely to have the same x-dependence. He then notes that if the 1 = 
configuration dominates as x — > 1, the observed ratio /F$ = 1/4 would follow. If we take 
the modern view that this high x behaviour will be controlled by the valence quarks, and 
note the quark model correlation between isospin and spin in the valence quark sector, this 
argument would also naively lead to the conclusion that u v ] (x) will dominate as x — ► 1. 
While correct, since Feynman's argument relies on the "wee" partons being uncorrelated 
with the leading quark, and so does not take into account the required antisymmetrization 
between the leading u quark and the "wee" u quark, its predictions for the full valence region 
are unclear. 



A more complete quark model argument is given in the papers of Close [|1J and Carlitz 
and Kaur ||14|| . They argued that 5'f/(6)-breaking changes Eq. (5) into 

p T= uudC A [^ t x° ud r + T xZ - t/f I x^ 1 ] (28) 



where X ° ud = y/l(U ~ It), xll =TT, and X Z = V§(U + IT) are the S = and two S = 1 
ud spin wavefunctions. For ip° = ip 1 = ip s , this wavefunction collapses to (5). Referring 
to hyperfine forces as driving the physics (which is equivalent to Feynman's assumption in 
this case), these papers posit that S'f/(6)-breaking leads to ip° ^ tp 1 , which would in turn 
lead to the relations u v j (x) = \v\(x) + |t> (x), u v [ (x) = ^vi(x), d v f (x) = |fi(x), 
and d v j (x) = \v\{x) in terms of distribution functions Vq and V\ associated with ip° and 
ip 1 , respectively, with V\/vq — > as i — > 1. This model thus also leads to F^/Ff — > j 
as x ^ 1 and it predicts that both A\{x) and A™( 1. While assuming 

ip° ^ if) 1 is very natural, since these diquark spin states are eigenstates of the hyperfine 
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interaction, this assumption is not consistent with the Pauli principle unless ip° and ip l 
have very special properties under the permutation group S3 or the wavefunction (28) is 
antisymmetrized. Thus, as with Feynman's argument, it is unclear what these models 

is 



predict for the full valence region. The closely related model of Close and Thomas |T5 
based on examining the energy of the spectator diquark after the deep inelastic scattering in 
a rest frame calculation of deep inelastic structure functions. Although, as pointed out by 
the authors, their calculation suffers from the fact that the diquark is a colored object which 
cannot have a well-defined energy, this calculation emphasizes the same physics and reaches 
the same conclusions as Refs. Given that the impact of the hyperfine interaction is 

implemented somewhat intuitively in this work, it is once again unclear whether the results 
presented are reliable for anything other than the x — > 1 behaviour. 

Although they do not use the hyperfine-perturbed quark model, the formalism required 
to deal explicitly with the fully antisymmetrized nucleon wavefunction seems to have first 



been applied to the valence quark spin distribution functions by Le Yaouanc et al. |L6 
They introduce an SU(6) intraband mixing between the ground state [56, + ] and the mixed 
symmetry [70, + ] in an attempt to account for the observed behaviour i^/Ff — > \ as 
x — > 1, i.e., d{x)/u(x) — > as x — > 1. This is precisely the kind of mixing introduced in (18) 
as required by color hyperfine interactions. (In fact, using this formalism makes calculations 
much simpler than in the uds basis, though perhaps less physically transparent.) They then 
make a prescription to boost this mixed wavefunction into the infinite momentum frame, fit 
the mixing angle to the data, and discuss the implications of such mixing to a wide range 
of phenomena. 

More recently, a number of authors |T7^p0| have addressed the connection between the 
hyperfine-perturbed quark model (either potential-based or bag-like) and the quark distri- 
bution functions. Most of these papers find the same two key effects I have emphasized here: 
axial current quenching by internal quark motion and u quark dominance as x — ► 1. 

Despite this extensive body of work [P~2"|-f2T)|] , it does not seem to be widely appreciated 
that the hyperfine-perturbed valence quark model makes quite clear predictions for the 
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asymmetries A^(x) and A r {(x) in the valence region. I attribute this state of affairs to 
the fact that this work has been very ambitious: most authors have attempted "absolute" 
calculations of structure functions. In doing so they encountered many obstacles, which 
forced them to a variety of assumptions, approximations, and "procedures" . The result is a 
wide range of predictions for the structure functions with apparent agreement only on their 
qualitative features. 



V. CONCLUSIONS 



In this paper I have shown that once it is assumed that the hyperfine perturbations of 
the quark model are responsible for the SU (6)-breaking observed in the structure functions, 
a very narrow band of predictions follows. In a broader context, I have argued that the 
extensive measurements and theoretical studies engendered by the "spin crisis" should be 
anchored in knowledge of whether the valence quark spin distributions are in fact anomalous. 
Thus whether the distributions described here prove to be correct when confronted with the 
data will be interesting, but not as important as the fact that such data will indicate whether 
the valence spin structure functions are in fact anomalous, and thus guide the search for 
where the resolution of the "spin crisis" is to be found. 
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